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SECTION  » 


INTRODUCTION 

Tho  recent  Interest  In  the  application  of  laminated  composite  materials 
to  forward  swept  wing  structural  design  has  led  to  studies  Indicating  that 
static  a>.me1ast1c  performance  may  be  enhanced  considerably  by  aeroelastic 
tallorii).;^  of  forward  swept  wing  structures.  Reference  1 describes.  In 
detail,  the  background  of  the  forward  swept  wing  concept  and, In  addition, 
presents  an  algebraic  expression  relating  the  divergence  speed  of  swept 
wings  to  the  geometric,  aerodynamic  and  structural  parameters  of  such  wings. 
The  present  report  describes  the  fundamental  behavior  of  swept  wings  con- 
stY,.:t<?.  with  composites  In  the  sub-critical  flight  regime. 

T spanwise  Ilf'  distribution  acting  upon  a wing  dictates  the  structural 
design  of  the  wing  support  structure.  Aeroelastic  effects  occur  because  the 
aerodynamic  loads  cause  deformation  of  a flexible  structure.  These  defor- 
mations cause  changes  in  the  shape  of  the  lifting  surface  that,  in  turn, 
lead  to  changes  In  the  applied  loads.  These  interrelated  effects  occur  until 
a deformed  equilibrium  state  is  reached  for  which  the  internal  forces  due 
to  deformation  balance  the  modified  aerodynamic  loads.  The  flexibility  of 
a lifting  surface,  because  of  its  effect  on  the  magnitude  and  distribution 
of  the  external  loads,  causes  changes  In  the  wing  lift-curve  slope,  wing 
bending  and  torsional  moments,  the  spanwise  center  of  nressure  and,  for 
swept  wings,  changes  in  the  aircraft  longitudinal  center  of  pressure. 
Aeroelastic  'Jeformatlon  and  its  effect  upon  spanwise  loading  is  recognized 
to  be  important  in  a number  of  situations  such  as:  wings  operating  at 
high  dynamic  pressures;  thin  wings;  swept  wings;  and  wings  designed  for 
low  wing  loading  [2]. 
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For  a given  wing,  the  effect  of  static  aeroelasticity  will  Increase 
with  flight  dynamic  pressure  since  the  lift  per  unit  change  In  angle  of 
attack  Is  proportional  to  dynamic  pressure.  It  Is  theoretically  possible 
to  encounter  a flight  condition  for  which  the  change  In  the  external  lift 
load  due  to  structural  deformation  exceeds  the  structural  restoring  forces 
within  the  structure  This  condition  produces  the  structural  Instability 
referred  to  as  wing  aemelas^lic  divergence. 

Modem  day  computeri;  .-d  computational  methods  now  exist  to  analyze 
the  effects  of  sta^^ic  aeroelasticity  upon  modern  con^slte  wings.  However, 
since  these  methods  r '<.sre  a reasonable  amount  of  knowledge  of  the 
structural  characteristics  of  the  wing  Itself,  they  are  fairly  cumbersome 
and  Inefficient  to  use  for  parameter  studies  whose  Intent  Is  to  develop 
Information  about  the  qualitative  oehavlor  of  laminated  comoosite  structures. 
A need  exists  to  determine  some  of  the  more  Important  consequences  and 
benefits  of  the  use  of  aeroelastic  tailoring  on  swept  wing  stmctures. 

The  present  study  extends  the  studies  presented  In  Reference  1 to  d1s> 
cuss  the  potential  for  utilization  of  laminated  composite  materials  to  en- 
hance the  pre-divergence  or  sub-critical  static  aeroelastic  characteristics 
of  swept  wings  In  general  and  forward  swept  wings  In  partlcul-r.  Specifi- 
cally to  be  considered  In  the  present  mport  are:  the  effects  of  laminate 
parameters  such  as  fiber  orientation  and  laminate  lay-up  on  divergence; 
load  redistribution; and  lateral  control  effectiveness  of  swept  wings. 

To  accomplish  these  objectives,  two  theoretical  models  were  formulated 
for  the  present  study.  The  first  of  these  models  uses  a laminated  box- 
beam  Idealization,  together  with  simple  strip  theory  air  loads,  to  study 


the  spanwise  lift  redistribution  problem  that  arises  for  flexible,  uniform 
planform,  swept  wings.  The  other  model  uses  a more  general  laminated  box- 
beam  idealization  of  the  wing  structure,  and  employs  a more  accurate  and 
versatile  representation  of  the  aerodynamic  loading  to  study  more  realistic 
configurations. 
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SECTION  II 


SUBCRITICAL  STATIC  AEROELASTIC  RESPONSE  - 
UNIFORM  PLANFORN  WING,  STRIP  THEORY  AIRLOADS 

The  differential  equations  of  equilibrium  for  a slender  wing  with 
moderate-to-high  aspect-ratio  may  be  formulated  easily  when  the  spanwise 
centerline  of  the  wing  box  is  restricted  to  be  a straight  line  and  the  de- 
formations of  the  wing  are  described  in  terms  of  elementary  Euler-Bernoulli 
beam  bending  theory.  This  approach  is  used  by  Diederich  and  Foss  [2]  to 
study  the  static  aeroelastic  behavior  of  metallic  swept  wings.  In  this  sec- 
tion this  elementary  approach  is  used  to  study  the  static  aeroelastic  defor- 
mation problem  for  composite  swept  wings.  The  resulting  equations  are  solved 
for  the  simple  case  in  which  the  geometrical  and  structural  properties  are 
constant  along  the  wing  span. 

1.  SLENOER  COMPOSITE  SWEPT  WING  EQUILIBRIUM  EQUATIONS 

The  differential  equations  of  equilibrium,  in  terms  of  bending  and 
torsional  deformations,  of  a slender  swef-t  wing  may  be  formulated  oy  con- 
sidering the  aerodynamic  loads  to  act  upon  chordwise  segments  of  the  wing, 
perpendicular  to  the  wing  reference  axis,  as  shown  in  Figure  1.  The 
structural  model  for  the  wing  assumes  the  loads  to  be  carried  by  a box 
beam  arrangement  constructed  of  laminated  composite  materials  forming  the 
upper  and  lower  face-sheets  of  the  wing  box.  The  reference  axis,  labelled 
as  the  y axis  in  Figure  1,  is  located  equidistant  between  the  front  and 
rear  edges  of  the  box  beam  and  lies  in  the  geometric  midplane  of  the  box, 
halfway  between  the  upper  and  lower  surfaces  of  the  cover  sheets. 

Reference  1 describes,  in  detail,  the  theoretical  development  of  the 
lamina ted- composite,  box-beam, structural  model.  The  governing  equations 


cf  equIllbrIuM  for  such  a wing  are  found  to  be  (note  that  ( )'  « d(  )/dy 
In  what  follows): 

(Elh-  - Ko')“  - p(y)  (1) 

(-Kh“  + 6Jo')*  * -t(y)  (2) 


Figure  1 - Geoinetry  of  idealized  swept  wing,  showing  a chordwise 
section  used  in  the  computation  of  aerodynamic  loads. 

The  structural  parameters  El,  GJ  and  K are  the  bearu  bending  stiffness, 
torsional  stiffness  and  bending-torsion  coupling  parameter,  respectively, 
and  are  derived  and  defined  in  Reference  1.  These  quantities  are  given  in 
Appendix  0.  The  variables  h(y)  and  a(y)  represent  the  upward  bending 
displacement  of  the  wing  reference  axis  and  the  nose-up  twisting  rotation 
of  the  wing  sections,  respectively.  The  term  p(y)  is  the  upward  load  per 
unit  length,  measured  along  the  y-axis,  while  t(y)  is  the  nose-up  torque 
per  unit  length,  also  measured  along  the  y-axis. 

Reference  3 describes  one  method  that  may  be  used  to  compute  the 
wing  loads  p(y)  and  t(y)  for  this  type  of  slender  wing  model.  For  this 
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calculation,  aerodynamic  strip  theory  Is  used  to  develop  the  following 
expressions: 


p(y)  = cc|f  q cos^A  + qca  cos^A(o  - rt&nA)  - nw  (3) 

X w 

r ? 2 2 

t(y)  = qcec^  cos*A  + ^ 

2 

+ qceaQ(a  - rtanA)  cos  A 


w% 

In  equations  (3)  and  (4),  the  term  c^  corresponds  to  the  two-dimension- 
al lift  coefficient  for  wing  sections  perpendicular  to  the  reference  axis 
for  a wing  with  no  flexibility, while  a^  refers  to  the  two-dimensional 
lift-curve  slope  of  the  same  section,  r Is  the  slope  of  the  bending  de- 
formation, defined  as  r = dh/dy,wh11e  nw  represents  the  distributed 
Inertia  load  per  unit  length  for  a wing  of  weight  w pounds  per  unit  length 
operating  at  a load  factor  n.  The  substitution  of  Equation  (3)  Into 
Equation  (l),and  Equation  (4)  Into  Equation  (2)  provides  the  differential 
equations  of  equilibrium  that  govern  the  static  aeroelastic  response  pro- 
blem. Equation  (1)  Is  the  bending  equation  and  can  be  nondimenslonallzed 
to  yield  the  following  relationship: 


3 2 

A^v  A^^  a„  cos  A 

O r . k + -I 2 (a  - rtanA) 


“T  ■ ~3 

dn  dn^ 


El 


= cos^A/EI  + nwt^/EI 


(5) 
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where  n " « - y/i  and  r * ^ • 

With  similar  nondimensional  definitions.  Equation  (2),  the  torsion  equation, 
becomes: 

— % - g + (a  - rtanA) 

dn"^  dn"^  GJ 


^ -qcei^  c]^  cos^A 
GJ 


(i  ♦ (|)(^)) 


(6) 


+ nwdt^/GJ 


Equations  (5)  and  (6)  are  coupled  together  by  the  dependence  of  the 
aerodynamic  load  upon  the  term  (a-rtanA)  and  by  the  structural  coupling 
parameter  K and  its  nondimensional  counterparts  k - K/EI  and  g = K/GJ. 

These  two  equations  can  be  combined  to  form  a single  equation  in  terms  of 
a new  variable  a^,  defined  as  = a - rtanA.  The  procedure  used  to  achieve 
this  combination  is  as  follcws: 

1)  Differentiate  the  torsion  equation.  Equation  (6)  with  respect 
to  n;  then,  multiply  the  result  by  the  factor  (1  - k tanA). 

2)  Multiply  the  bending  equation.  Equation  (5),  by  the  factor 
(g  - tanA). 

3)  Add  the  two  above  results  together  to  obtain  the  following 
equation; 

f. 


a'r  + aal  - ba.  = f^’ 


e e "“e  l 2 
with  the  definitions  of  terms  a,b,f^f2  given  as 

2 , 2 

/qcer 
a 


(7) 


^ |qce£  a^cosAj  ^ „ tanA  j 


(8) 

19) 


7 


The  boundary  conditions.  In  terms  of  the  variable  o^,  are  found  by 
first  noting  that  both  r (the  bending  slope)  and  a are  zero  at  the  effec- 
tive root  of  the  wing  (n  = 1).  This  leads  to  the  boundary  condition 

ag(r,  = l)  = 0 (12) 

The  fact  that  the  bending  moment  and  the  twist  are  zero  at  the  effective 
wing  tip  leads  to  the  equation 

cte'(n  = 0)  = 0 (13) 

The  wing  tip  zero  shear  condition,  when  combined  with  the  torsion  equation, 
evaluated  at  n = 0,  leads  to  the  following  result: 

ag"(0)  + a ag(0)  = f^(0)  (14) 

2.  AEROELASTIC  EFFECT  EQUATION  SOLUTION 


To  solve  the  wing  deformation  problem  defined  by  Equation  (7)  and 
its  associated  boundary  conditions,  we  restrict  ourselves  to  an  untwisted, 
uniform  planform  wing  whose  pitch  attitude  is  specified.  In  this  case, 
f^(n)  and  f2(n)  are  constants.  It  is  shown  in  References  2 and  3 that  the 
solution  to  this  problem  is  given  by 


The  function  fj(n)  is  defined  as 


f,(n) 

’ w 


(15) 
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fo(n) 


4B 

x= J 

9 


( A 6^  . A 

[l?7T7/ 


cosYn  + 


3 B^  - 

— jr— 

9 B^  Y 


sinYn 


The  terms  B and  y that  appear  in  Equation  (16)  are  obtained  from  the 
knowledge  that  -2B  and  6 + 1y  are  roots  of  the  following  characteristic 
equation,  obtained  from  Equation  (7), 

r^  + ar  - b = 0 (17) 

If  the  parameters  c-j  and  C2  are  defined  as 


Cl  = Mr  + 


2 ■"Vr^?7 


r ^ JT 


then  B and  y may  be  written  as 


e=-|(Ci+C2)  (20) 

Y = (c^  - Cg)  (21) 

Equation  (15)  provides  Information  necessary  to  compute  the  ratio 
between  the  total  lift  developed  by  a flexible  composite  wing  and  the  lift 
developed  by  a similar,  but  Inflexible  wing.  Under  the  strip  theory  re- 
striction, the  lift  of  the  rigid  wing  is  found  to  be 

“ 9cc^  ^ cos^A  - nwt  (22) 

The  lift  on  the  flexible  wing  Is 


L = L. 


(qc  a^  cos^A)  Og  dy 
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Performing  the  necessary  integration  in  Equation  (23)  and  dividing  by 


yields  an  expression  for  the  flexible-to-rlgid  wing  lift  ratio: 
k_  -2by6~^^  + 2BYe^cosY  (3P^  y^)e^  siny 


r 4B^e”^®  + e^[(5B^  + y^)  cosy  + (3B^  - By^)  siny] 


(24) 


Equation  (24)  is  an  expression  for  the  lift  effectiveness  of  the  wing. 
Note  that  when  the  denomitiator  becomes  zero,  then  the  wing  will  diverge. 

The  ratio  of  the  bending  moment  at  the  wing  root  for  the  flexible 
wing  to  that  computed  for  the  rigid  wing  is  found  to  be 


!L=  2 

«r 


ye'^^  + e^(-  ycosy  + 3g  siny) 

4B^ye~^^  + e®[(5B^Y  + y^)cosy  + (3B^  - BY^)siny] 


(25) 


The  bending  moment  is  measured  about  the  effective  root  of  the  wing  as  shown 
in  Figure  1. 

The  twisting  moment  ratio  expression  is  identical  to  that  calculated 
for  the  ratio  L/L^.  The  spanwise  center  of  pressure,  Y^p,  measured  along 
the  reference  axis,  can  be  obtained  by  combining  Equations  (24)  and  (25) 
and  reads  as  follows: 

’'cp4(^](r)  <“> 

3.  DISCUSSION  OF  RESULTS 

a.  The  Change  in  the  Position  of  the  Center  of  Pressure 

The  change  in  the  spanwise  position  of  the  center  of  pressure, due  to 
wing  flexibility, affects  both  the  overall  directional  stability  of  the 
aircraft  and  the  wing  root  stresses.  An  outboard  shift  of  the  center  of 
pressure  results  in  an  increase  in  stress  levels  near  the  wing  root.  The 
equations  obtained  previously  may  be  used  to  compute  the  change  in  the 
spanwise  location  of  the  center  of  pressure  position  along  the  swept 
reference  axis.  Reference  2 presents  a study  of  this  change  for  uniform 
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planform  metallic  wings,  as  a function  of  two  parameters.  These  parameters 
are  related  to  dynamic  pressure  and  to  the  sweep  angle.  The  present  study 
presents  similar  results  by  suitably  redefining  parameters  to  conform  to 
the  composite  wing  analysis.  The  two  parameters  are  the  parameter  a,  defined 
In  Equation  (8),  and  the  ratio  of  b/a  defined  as  follows: 


The  parameter  a Is  referred  to  as  the  dynamic  pressure  parameter,  while 
b/a  Is  termed  the  sweep  parameter.  Figure  2 shows  the  Influence  of  these 
parameters  upon  the  movement  of  the  center  of  pressure  (CP)  along  the  swept 
reference  axis.  The  parameter  ay*  Is  equal  toAYj.p/)i/2  , and  positive  ay* 
corresponds  to  an  outward  CP  movement  while,  conversely.  Inboard  CP  movement 
Is  denoted  by  negative  values  of  ay*. 

Labelling  the  upper  right  quadrant  of  Figure  2 as  Quadrant  1,  and 
numbering  the  other  quadrants  consecutively,  counterclockwise,  we  see  that 
values  of  a and  b/a  that  lie  In  Quadrant  1 correspond  to  wings  for  which 
divergence  Is  very  unlikely  to  occur.  In  addition,  the  center  of  pressure 
movement  Is  Inboard  for  values  of  b/a  greater  than  ...out  4. 

In  Quadrant  2,  divergence  Is  very  likely,  except  for  wings  falling 
Into  a very  narrow  range  of  a and  b/a  values.  In  addition,  shifts  In  the 
CP  location  are  outboard  so  that  unfavorable  effects  of  aeroelasticity 
upon  wing  stresses  may  be  expected  In  this  region. 

In  Quadrant  3,  divergence  Is  Impossible,  while.  In  Quadrant  4,  diver- 
gence may  be  possible  for  certain  c mbinatlons  of  a and  b/a. 
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DYNAMIC  PRESSURE  PARAMETER, 


3 


SWEEP  PARAMETER,  b/a 

Figure  2-  Static  aeroelastic  characteristics  of 
a uniform  property  swept  wing,  showing 
spanwise  center  of  pressure  change  Ay* 
as  a function  of  the  parameters  a and 
b/a 
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Reference  2 notes  that  subsonic*  metallic  sweptback  wings  (k  * g « 0) 
are  described  by  values  of  a and  b/a  lying  In  the  first  quadrant,  while 
subsonic,  metallic  sweptforward  wings  fall  Into  Quadrant  2.  For  a to  be 
negative  for  a metallic  wing,  the  chordwise  aerodynamic  center  must  lie 
behind  the  wing  elastic  axis.  This  situation  Is  usually  associated  with 
supersonic  flight.  Quadrant  3 Illustrates  sweptforward  metallic  wing  be- 
havior in  the  supersonic  flight  regime,  while  Quadrant  4 displays  metallic 
sweptback  wing  behavior  at  supersonic  speeds. 

When  the  wing  Is  constructed  of  conq)os1te  materials,  the  likely  effect 
of  sweep  on  aeroelastic  behavior  Is  not  as  clear  as  It  Is  for  metallic  wings. 
The  elastic  coupling  between  bending  and  torsion  Introduces  a new  parameter 
into  the  definitions  of  a and  b/a.  For  a particular  value  of  wing  sweep 
angle  and  dynamic  pressure,  the  values  of  a and  b/a  may  be  modified  by 
tailoring  of  the  composite  laminate. 

b.  Example  Case 

Figures  3 and  4 are  presented  to  illustrate  the  potential  effect  of 
tailoring  laminate  fibers  for  aeroelastic  performance.  Consider  first 
Figure  3,  in  which  the  parameter  a is  plotted  versus  the  sweep  parameter 
b/a  for  an  example  wing  with  zero  sweep.  The  wing  is  one  of  constant  chord 
and  uniform  elastic  properties  along  the  span.  The  laminated  wirg 
structure  itself  has  a laminate  thickness- to-box  beam  depth  ratio  of  1:10. 
The  laminate  cover  sheets  are  composed  of  graphite-epoxy  with  lOX  of  the 
laminate  fibers  oriented  along  the  x-axis  (see  Figure  1);  25S  of  the  fibers 
are  along  the  + 45®  direction  (with  respect  to  the  y-axis)  and  65X  of  the 
fibers  are  oriented  at  a variable  angle  e,  measured  with  respect  to  the 
x-axis. 
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The  dynamic  pressure  value  chosen  corresponds  to  q equal 
to  the  divergence  q of  an  unswept  wing  with  a laminate  for  which  e 90°. 
The  plot  of  a versus  b/a  Is  symmetrical  about  the  vertical  axis  and  ex* 
hibits  a cusp.  At  this  point,  laminates  with  e = 0°  and  e - 90°  exhibit 
the  same  static  aeroelastic  behavior.  Referring  back  to  Figure  2,  It 
Is  seen  that  the  effect  of  sweeping  the  e fibers  forward  Is  to  keep  the 
wing  static  aeroelastic  characteristics  In  Quadrant  1.  The  favorable 
performance  of  the  wing  structure  In  this  region  Is  readily  apparent. 

Figure  4 shows  a plot  of  a versus  b/a  for  two  uniform  property 
wings,  one  with  15°  of  forward  sweep  and  the  other  with  15°  of  rearward 
sweep.  The  two  plots  are  seen  to  be  symmetrical  about  the  line  b/a  equal 
to  zero.  For  the  A = -15°  wing,  a cusp  point  Is  located  to  the  left  of 
the  origin  and  coincides  with  a wing  with  0 fibers  at  0°  and  approximately 
73°.  The  A = +15°  wing  has  a similar  cusp  when  e *0°  and  e =107°.  For 
this  wing.  It  Is  still  possible  to  design  the  laminate  such  that  diver- 
gence will  not  occur,  even  for  15°  of  forward  sweep.  Sweeping  the  wing 
beyond  15°  of  forward  sweep  Is  found  to  move  the  cusp  further  to  the 
lert  and  to  further  disort  the  figure. 
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SECTION  III 

MODERATE  TO  HIGH-ASPECT-RATIO  WING  OF  ARBITRARY  PLANFORM  IN  SUBSONIC  FLOW 
This  section  discusses  a computation  method  to  determine  the  sub- 
critical  static  deformation  of  a flexible,  moderate- to-high-aspect-ratio 
wing  in  a subsonic  flow.  The  wing  planform  itself  may  have  an  arbitrary 
shape,  while  the  structural  model  is  idealized  to  be  a nonuniform  box 
beam  constructed  of  composite  material.  The  wing  aerodynamic  forces  are 
computed  by  use  of  a modified  Weissirger  L-method,  valid  ct  subcritical 
Mach  numbers.  The  basic  computational  schemp  is  based  upon  that  found  in 
NACA  TN  3030  [4],  with  the  exception  that  the  method  has  been  updated  to 
include  composite  material  coupling  effects  such  as  those  discussed  pre- 
viously in  Section  II  of  this  report. 

1.  THE  AEROELASTIC  EQUATIONS  IN  MATRIX  FORM 

The  subcritical  static  aeroelastic  analysis  problem  can  be  reduced  to 
finding  the  solution  to  a set  of  simultaneous  linear  algebraic  equations. 
The  first  step  in  such  an  analysis  is  the  formulation  of  a relationship 
between  the  airload  distribution  along  the  wing  (lift  per  unit  length)  and 
the  streamwise  angle  of  attack  distribution.  The  formulation  of  the 
problem  involves  the  idealization  of  the  continuous  wing  into  a finite 


E 

\ 


number  of  discrete  panels.  Figure  5 illustrates  a typical  wing  idealization. 
The  lift  per  unit  length  p^  is  computed  at  the  center-span  of  the  panels, 
while  the  streamwise  angle  of  attack  a-  (positive  nose-up)  is  measured 
at  the  same  points.  The  lift  acting  on  each  panel  is  then  equal  to 

Pih.. 

Denoting  the  lift  intensity  at  a finite  number  of  positions  or  panels 
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(a)  ACTUAL  PLA^  -ORM 


Fiqure  5-  Wing  aerodynamic  model  idealization; 

(a)  actual  planform,  (b)  planform 
idealization 


along  the  wing  as  a vector  {p.}  and  the  streamwise  angle  of  attack  dis- 
tribution at  the  same  positions  as  a linear  relationship  can  be  writ- 
ten to  relate  {p^},  {a.j}and  the  dynamic  pressure,  q.  Expressed  in  standard 
matrix  not'tion,  this  equation  is 

l/n[A^j]{Pj}  = (28) 

The  square  matrix  [A. J is  the  matrix  of  aerodynamic  influence  coefficients; 

® J 

the  computation  of  this  matrix  Is  discussed  in  Appendix  a.  The  distribution  of 
angle  of  attack  {a^-}  along  the  wing  is  composed  of  three  types  of  com- 
ponents, and  is  expressed  as  follows: 

= {a^}  + {ag}  + {a^}  (29) 

The  portion  of  {o}  denoted  as  {a^}  is  a rigid  body  angle  of  attack  and  oc- 
curs either  because  of  the  incidence  angle  of  the  entire  aircraft  or  be- 
cause of  a predetermined  initial  or  "built-in"  geometric  angle  of  attack 
distribution  along  the  wing.  The  vector  {og}  represents  the  angle  of  at- 
tack distribution  caused  by  the  deformation  of  the  flexible  wing  structure. 

The  {a^}  distribution  represents  that  part  of  the  angle  of  attack  distri- 
bution due  to  apparent  or  actual  aerodynamic  twists,  such  as  those  caused 
by  aileron  deflection  or  steady-state  angular  velocity  due  to  aircraft 
roll  motion. 

The  angle  of  attack  distribution  due  to  structural  deformation  is 
computed  from  the  equation 

{as>  = [C^j]{Pj}  (30) 

where  C..  represents  an  element  of  the  wing  flexibility  matrix.  For  the 

• J 

computation  of  this  matrix  it  is  assumed  that  the  wing  semi -span  is 
clamped  at  its  root  and  free  at  the  tip.  The  details  of  the  formulation  of 


this  matrix  for  a wing  with  a composite  material  wing  box,  idealized  as 
a beam,  are  discussed  in  detail  in  Appendix  A.  These  coefficients  had  not 
been  derived  previously  in  the  open  literature. 

By  combining  Equations  (28),  (29)  and  (30),  a basic  relationship,  re- 
lating loads  to  angles  of  attack,  is  found  to  be  as  follows: 

CBijKpj}  = {Or}  + {ac>  (31) 

where 

CB^j]  = [l/q[A^j]  - [C^j]]  (32) 

Computational  efficiency  is  achieved  when  the  general  wing  airload 

distribution  is  subdivided  into  two  basic  types:  those  loads  due  to 
maneuvers  resulting  in  symmetrical  lift  distributions;  and  those  maneuvers, 
such  as  may  result  from  aileron  actuation,  resulting  in  asymmetrical  lift 
distribution.  When  this  subdivision  of  load  sets  is  used,  only  one  wing 
semi-span  need  be  studied.  In  Reference  4,  and  in  this  report,  the  left 
wing  is  used.  Thus,  the  matrix  [B^j]  is  computed  for  half  the  wing  span, 
once  for  the  symmetrical  maneuver  (e.g.,  the  1-g  load)  and  once  for  the 
asymmetrical  maneuver  (e.g.,  the  aileron  effectiveness  problem). 

The  matrix  [B-.],  as  defined  in  Equation  (32),  is  a function  of  the 

' J 

flight  dynamic  pressure,  q.  The  solution  for  the  airload  distribution 
{p.}  requires  computation  of  the  inverse  of  Wing  configurations, 

for  which  [B^j]  is  singular  at  certain  critical  values  of  q may  exist. 

These  critical  values  of  q,  denoted  as  q^,  correspond  to  the  clamped  wing 
static  divergence  dynamic  pressures.  For  a symmetrical  maneuver,  the  load 
distribution  along  the  wing,  is  given  by  the  expression 

{p.}  = [B. 


(33) 


(34) 


The  total  lift,  L,  of  the  Ming  is 


1.  = 2 LIJ  {Pi>  * 2 [IJ 


If  {a^}  Is  given  by  the  equation 

{op}*  CqCI) 

then  Equation  (34)  is  written  as 

L = 2aQLlJ[B.jr‘tl} 

The  lift  curve  slope  of  the  wing  is  written  as 

C^^  = L^q  S ) 


= (2/qS)Ll  J[B.jr^{l} 


where  S is  defined  to  be  the  total  wing  pi  anf  ora  area.  A two-dimensional  sec- 
tional lift-curve  slope  is  included  in  the  input  parameters  that  must  be 
known  before  the  coefficients  can  be  calculated.  The  effect  of  Mach 
number  on  the  wing  lift-curve  slope  may  be  accounted  for  by  modification  of 
the  2-D  sectional  lift  curve  slopes  for  each  section.  The  reader  is  re- 
ferred to  Reference  4 for  a complete  discussion  of  this  modification. 

The  effectiveness  of  aileron  deflection  in  generating  roll  may  be  com- 
puted by  a method  similar  to  that  outlined  for  the  symmetrical  airload 
computation.  Since  ailerons  are  deflected  asymmetrically  to  initiate  and 
sustain  rolling  motion,  the  lift  distribution  that  results  from  roll  is  asym- 
metrical. The  asymmetrical  Influence  of  one  wing  semi -span  upon  the  other  is 
accounted  for  by  proper  computation  of  the  aerodynamic  influence  coefficients. 
This  matrix  of  coefficients  for  the  asymmetrical  lift  case  is  denoted  as 


The  equation  of  equilibrium  for  an  aircraft  undergoing  a steady>state 


roll  may  be  written  as: 

1 (Pj)  • [C,j]  (Pj)  * q[E,j]  (« 

^ ,3o,  + Pbl^^l 
^ ^ 5T|“) 


(39) 


Referring  to  the  right-hand  side  of  Equation  (39),  the  first  term  corres- 
ponds to  the  angle  of  attack  distribution  due  to  structural  flexibility. 

The  second  tenn,  involving  the  square  matrix  appears  because  the 

deflection  of  an  aileron  causes  a change  in  the  pitching  moment  of  the 
various  wing  sections.  This  change  in  pitching  moment  produces  a change  in 
the  angle  of  attack  due  to  structural  deformation  and  is  proportional  to  q. 
The  detailed  development  of  the  matrix  is  discussed  in  Appendix  C. 

The  third  term  in  Equation  (39)  accounts  for  the  apparent  change  in 
angle  of  attack  of  a wing  section  that  occurs  with  aileron  deflection. 
Finally,  the  term  proportional  to  Pb/2V  represents  the  damping- in-roll  term. 
This  latter  term  appears  because  rolling  motion  of  the  wing  (positive  if 
the  left  wing  moves  down)  produces  an  apparent  change  in  angle  of  attack 
at  a position  y.  on  the  left  wing,  given  by 

Ao^  = Py/V  (40) 

When  the  aircraft  is  in  a steady-state  roll  condition,  the  equation  for 
moment  equilibrium  becomes: 

Ly^J{p^h.}  = 0 (41) 

Equation  (41)  is  used  to  solve  for  the  amount  of  roll  rate  caused  by 
an  aileron  deflection  {6^}  along  the  left  wing,  with  a corresponding  de- 
flection {-6^-}  along  the  right  wing. 


Let 


= Wi/*’]  - t'lj]  («) 

while  {6^}  = 6^  {4,}  (43) 

Equation  (43)  allows  segmented  ailerons,  although,  in  the  usual  case,  this 
vector  would  be  composed  of  values  of  0 and  -1,  corresponding  to  no  aileron 
or  an  "up  aileron",  respectively,  at  a wing  station. 

Using  the  definitions  in  Equations  (41),  (42)  and  (43),  the  following 
relationship  is  found 

(Pb/2V)/«„  = — (44) 

’3  <-r> 

The  quantity  (Pb/2V)/6^  is  termed  the  aileron  effectiveness.  The  aileron 
effectiveness  is  a measure  of  the  effectiveness  of  lateral  control  devices 
and  represents  the  wing-tip  helix  angle  per  radian  of  asymmetrical  aileron 
displacement.  When  the  aileron  effectiveness  is  a positive  quantity, 
positive  roll  control  is  possible.  Since  aileron  effectiveness  is  a func- 
tion of  q,  there  may  be  situations  for  which  q is  large  enough  to  cause  con- 
trol reversal,  that  is,  the  aileron  effectiveness  may  be  negative.  The 
aileron  reversal  speed  corresponds  to  the  speed  at  which  the  aileron  ef- 
fectiveness is  zero. 


2.  AN  EXAMPLE  CASE 


To  illustrate  the  use  of  the  computational  method  described  in  the 
previous  section,  an  example  wing  was  developed.  The  wing  planform  has 
a linear  taper  and  a tip-to-root  chord  ratio  of  0.50;  the  reference  axis 
is  swept  forward  30**.  The  chordwise  sections  of  the  wing  are  geometrically 
similar:  this  leads  to  a situation  in  which  the  bending  stiffness,  torsional 
stiffness  and  coupling  parameter  vary  along  the  reference  axis  as  follows. 

£I  = f^  EI^;  6J  = f^  6J^;  K = f^  (45) 

where  f = 1 - n (1-x) 
n = y/i 

X = taper  ratio  = 
and 

( )p=  root  value 

The  wing  structure  itself  consists  of  thin,  laminated  cover  sheets, 
t p and  bottom.  These  laminates  are  composed  of  graphite-epoxy  lamina  and 
have  a layup  identical  to  those  described  for  the  example  in  Section  2.3(b). 
The  wing  planform  has  a reference  chord  of  45  inches  and  a semi -span,  b/2, 
of  120  inches. 

Figure  6 shows  a plot  of  nondimensional  divergence  velocity  at  sea 
level  in  incompressible  flow.  The  reference  velocity,  Vp,  corresponds  to 
the  divergence  speed  of  the  wing  when  the  variable,  or  e,  fibers  lie 
along  the  reference  axis,  90**.  The  maximum  divergence  velocity  is  seen 
to  occur  when  e is  approximately  110®,  or  20®  forward  of  the  reference  axis. 
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At  this  position,  the  wing  structure  is  stiffest  from  the  standpoint  of 
divergence. 

To  illustrate  typical  sub-critical  aeroelastic  characteristics  of  the 
laminated,  forward  swept  wing,  a situation  is  examined  for  which  the  dynamic 
pressure  is  20%  of  the  divergence  dynamic  pressure  of  the  wing  with  e fibers 
oriented  at  90°.  The  effect  of  fiber  angle  upon  aeroelastic  stiffness  can 
be  seen  in  Figure  7.  In  this  figure,  the  parameter  Ay*  is  plotted  against 
0.  The  parameter  ^y*  represents  the  spanwise  movement  (positive  outward) 
of  the  wing  center  of  lift  position,  from  the  center  of  lift  position  on  the 
rigid  wing,  divided  by  the  wing  semi -span.  Large  values  of  this  tiy* 
parameter  are  associated  with  large  positive  changes  in  the  root  bending 
moment.  This  is  true  because,  for  a given  flight  condition,  the  lift  re- 
nains  constant.  Near  the  e = 110°  position,  the  outward  movement  of  the 
center  of  lift  is  very  slight,  while,  at  e * 40°,  a relatively  large  out- 
ward shift  in  the  center  of  lift  can  be  seen. 

Turning  to  aileron  effectiveness.  Figure  8 illustrates  the  behavior 
of  the  flexible-to-rigid  ratio  of  Pb/2V/6Q  versus  fiber  angle  e for  the 
example  wing.  Two  cases  are  considered.  The  first  case  illustrates  the 
aileron  effectiveness  of  the  wing  when  the  ailerons  are  located  between 
0.30  < 2 y/jj  < 1.0,  while,  in  the  second  case,  the  ailerons  lie  in  the 
region  0.30  < 2y/^  < 0.80. 

For  Case  1,  the  maximum  aileron  effectiveness  is  achieved  when  e = 125°, 
while,  for  Case  2,  maximum  effectiveness  occurs  at  e = 115°.  These  e angles 
are  slightly  different  than  the  e = 110°  angle,  at  which  the  maximum  di- 
vergence speed  occurs. 
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forward  sweep,  as  a function  of  e.  Each  curve 
is  normalized  with  respect  to  the  effectiveness 
of  the  ailerons  on  a rigid  wing. 
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SECTION  IV 
CONCLUSIONS 

This  report  has  detailed  the  results  of  an  investigation  of  the  sub- 
critical  static  aeroelastic  characteristics  of  swept  wings.  In  addition, 
a matrix  method  of  analysis  is  detailed  which  will  provide  an  accurate,  but 
low-cost,  representation  of  the  aeroelastic  behavior  of  fibrous  composite 
swept  wings  in  subsonic  flow.  The  theories  developed  have  been  used  in  a 
limited  nun^er  of  cases  to  illustrate  the  potential  benefits  to  be  gained 
from  aeroelastic  tailoring  of  forward  swept  wings. 

Several  important  conclusions  may  be  drawn  from  this  analysis: 

(1)  The  lift  coefficient  of  a flexible  wing  can  be  modified  greatly  by 
the  judicious  choice  of  laminate  orientation.  Similarly,  the  wing 
spanwise  center  of  pressure  location  can  be  moved  significantly  by 
aeroelastic  tailoring. 

(2)  The  lateral  control  effectiveness  of  the  wing  may  be  altered  by  lamin- 
ate tailoring.  However,  maximum  aileron  effectiveness  does  not  occur 
when  the  divergence  speed  is  maximized. 

(3)  The  matrix  method  developed  has  the  potential  for  quickly  analyzing 
the  effects  of  laminate  construction  upon  the  aeroelastic  character- 
istics of  the  wing  at  low  cost. 

Finally,  the  techniques  developed  in  this  report  have  been  shown  to  pro- 
vide an  effective  means  of  further  exploring  the  influence  of  the  various 
parameters  upon  the  static  aeroelastic  characteristics  of  laminated,  swept 
wing  structures. 
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APPENDIX  A 

THE  DEVELOPMENT  OF  THE  WING  FLEXIBILITY  MATRIX 


The  wing  structural  Idealization  used  throughout  this  report  consists 
of  a box  beam  whose  flexibility  Is  entirely  due  to  the  presence  of  upper 
and  lower  cover-sheets  constructed  of  laminated  conq>os1te  material.  This 
appendix  details  the  development  of  the  flexibility  matrix  for  this  type 
of  wing  structural  model.  The  laminated  beam  theory  upon  which  the  for> 
mulation  is  based  is  discussed  in  detail  in  Reference  1. 

An  arbitrary  planform  wing  is  shown  in  Figure  A-1  to  illustrate  the 
geometry  used  in  the  flexibility  matrix  development.  The  wing  semi-span  is 
divided  into  a finite  number  of  streamwise  panels  of  arbitrary  width  h^.. 

The  continuously  varying  spanwise  airload  distribution  is  replaced  by  a set 
of  discrete  running  loads,  each  of  which  acts  upon  a streamwise  panel  and  is 
of  constant  magnitude  over  the  panel  on  which  it  acts. 

The  loads  acting  on  the  wing  panels  cause  wing  bending  and  torsion  de- 
formation. The  bending  and  torsion  moments  caused  by  the  applied  loads  are 
computed  in  local  coordinate  systems  located  on  each  panel  as  indicated  in 
Figure  A-2,  This  local  reference  system  is  swept  at  an  angle  to  the 

horizontal,  with  sweepback  being  positive.  The  elastic  properties  of  the 
panel  are  computed  with  respect  to  the  local  or  swept  axis. 

The  development  of  the  flexibility  matrix  for  the  composite  wing  is  based 
upon,  and  closely  parallels,  the  similar  development  of  a flexibility  matrix 
for  a metallic  swept  wing  detailed  in  Reference  4.  The  geometry  of  the 
structural  skeleton  of  the  wing  model  is  shown,  for  three  outboard  wing 
segments,  in  Figure  A-3.  Note  that  this  represents  the  left  wing  of  the 
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Figure  A-1  Four  panel  wing  model  showing: 

spanwise  panel  dimension  h^;  sweep 
angle  Ai;  and  segmented  aerodynamic 
centerlines  and  structural  reference 
axis. 
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Figure  A- 3 Plan  view  of  outboard  wing  sections. 

Geometrical  distances  are  positive  as 
shi 


aircraft.  The  definition  of  the  symbols  appearing  in  Figure  A-3  is  as 

follows: 

= total  lift  of  a spanwise  increment  of  wing  having  span  h^.  Panel 
numbering  is  sequential,  beginning  with  the  wing  tip. 

e^  = streamwise  distance  from  the  aerodynamic  center  of  the  panel  to  the 
structural  reference  point,  positive  when  the  reference  point  lies 
aft  of  the  aerodynamic  center. 

f.  = streamwise  distance  from  the  structural  reference  oointona  tarticular 
panel  to  the  reference  point  at  the  next  inboard  panel.  This  quantity 
is  positive  when  the  inboard  reference  point  lies  forward  of  the 
outboard  point. 

d^.  = lateral  distance  between  the  structural  reference  point  of  a panel 

and  the  structural  reference  point  of  the  next  inboard  panel. 

= sweep  angle  of  the  structural  reference  axis,  positive  for  sweep- 
back. 

= rolling  moment  at  a structural  reference  point  due  to  the  lift  of 
all  panel  segments  outboard  of  this  point,  positive  when  it  tends 
to  raise  the  left  wing  tip.  See  Figure  A-2. 

My.  = pitching  iroment  about  the  y-axis  at  a structural  reference  point, 
positive  nose-up.  See  Figure  A-2. 

M.  = beam  bending  moment,  at  a structural  reference  point,  about  an  axis 
perpendicular  to  the  local  axis^positive  when  the  upper  surface  is 
placed  in  compression. 

T.  = torsional  moment,  at  a structural  reference  point,  about  the  local 
axis,  oositive  when  it  tends  to  twist  the  leading  edge  up. 
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The  structural  reference  line  of  the  wing  is  composed  of  a set  of 
interconnected  straight  line  segments.  The  reference  axis  lies  midway  be- 
tween the  front  and  the  rear  edges  of  the  structural  box  beam.  Since  the 
load  intensity  is  assumed  constant  over  each  panel,  the  force  acting  on 
each  panel  at  the  reference  point  is 

L.  = p.h^  (A-1) 

where  p.  is  the  running  lift  per  unit  spanwise  length  (p.  = p^(y)).  As  a 
result  of  the  lift  loads  acting  on  each  panel  the  following  bending 
moments  and  torques  arise  [4], 


where 


{m}  = [[cosA]  [r^]  - [sinA]  [u]]  {l} 

|t|  = [CsinA]  [r^3  + [cosa]  [u]j 


(A-2) 

(A-3) 

(A-4) 


(A-5) 
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The  objective  of  the  flexibility  matrix  development  is  to  define  a 
relationship  between  the  load  intensities  p^  on  each  wing  panel  and  the 
streamwise  angle  of  attack  on  each  panel,  This  development  is  accom- 
plished through  the  energy  methods  and  the  use  of  the  dunn\y  unit  load 
method.  For  the  linear  elastic  laminated  beam  model  discussed  in  Reference 
1,  the  present  study  has  determined  that  the  strain  energy,  expressed  in 
terms  of  the  bending  moment  M and  twisting  moment  T,  is 

-i 

^ 1 ut  fitrru  Tt  I 

(A-12) 


U 


2j  1 - kg 


. 2K7K  . T^  I 

El  EIGJ  6J  j 


where  El,  GJ  and  K are  the  bending  stiffness,  torsional  stiffness  and 
coupling  parameter  values,  respectively  for  the  beam  and 

k = K/EI  g = K/GJ  (A-13a,b) 


ds  = incremental  distance  along 
the  reference  axis 


(A-13c) 


Note  that  all  values  are  computed  with  respect  to  the  structural  reference 
axis.  The  streamwise  angle  of  attack  distribution,  a(s),  along  the  wing, 
resulting  from  the  bending  moment  and  twisting  moment  distributions,  may 
be  found  by  using  Castigliano's  Theorem  to  derive  the  following  relationship 
for  a. 


<*-») 

The  terms  m and  t in  Equation  (A-14)  are,  by  definition,  the  bending  moment 
and  twisting  moment,  respectively,  that  arise  when  a unit  pitching  moment 
(positive  nose-up)  is  applied  at  a position  on  the  wing  where  a is  to  be 
computed.  Note  that  the  unit  pitching  moment  is  applied  perpendicular  to 
the  free  stream  while  m and  t are  computed  along  the  reference  axis.  Equation 
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(A-14)  may  be  used  to  compute  a in  matrix  form.  The  vector  |o|  represents 
the  set  of  streamwise  angles  of  attack  computed  at  the  spanwise  midpoint 
of  each  wing  panel.  The  values  of  El,  6J  and  K are  taken  to  be  constant 
for  each  panel  and  are  those  values  computed  at  the  panel  reference  point. 

The  general  procedure  for  deriving  +he  flexibility  matrix  for  an  iso- 
tropic metallic  wing  is  detailed  in  Reference  4.  Following  the  same  pro- 
cedure for  the  laminated  beam  model,  the  flexibility  matrix  [C^  •]  is  defined 
such  that  the  lift  distribution  |pj}is  related  to 

The  expression  for  [C..]  reads: 

■ J 


[C,j]  ’ [I„] 


0 

0 0 

h(g  - tan a) 
EI(1  - kg) 

[cosA][r^]  - [sinA][u] 

0 

0 0 

h(l  - ktanA) 

[sinAlCrg]  + [cosa][u] 

W(1  - kg) 

0 

[h] 

0 

[h] 


(A-16) 


where 


1 

2 


1 


0^11 
0 0 1 

0 0 0 |- 


(A-17) 


0 — ^ 

In  Equation  (A-16),  note  the  [ ] refers  to  a diagonal  matrix.  If 
k and  g are  zero,  the  expression  given  in  Equation  (A-16)  reduces  to  that 
derived  in  Reference  4. 


I 
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APPENDIX  B 

THE  AERODYNAMIC  INFLUENCE  COEFFICIENT  MATRIX 

This  appendix  reviews  the  theoretical  development  of  the  aerodynamic 
Influence  coefficient  matrix  [A^^]  used  In  this  study.  The  development  Is 
taken  from,  and  Is  Identical  to,  that  developed  In  Reference  4. 

The  spanwise  lift  distribution  on  an  aircraft  wing  varies  In  the  lateral 
direction  and  may  be  visualized  as  resulting  from  a system  of  horseshoe  vortices, 
each  of  constant  strength,  r^,  placed  along  the  wing  as  Indicated  In  Figure  B-1. 
The  double  arrows  Indicate  that  the  sense  of  the  circulation  around  each  vortex 
line  segment  Is  given  by  the  right  hand  rule. 

With  the  type  of  aerodynamic  model  shown  above,  the  net  strength  of  the 
trailing  vortex  at  any  point  on  the  wing  span  Is  numerically  equal  to  the  rate 
of  change  of  strength  of  the  bound  vortex  In  the  spanwise  direction.  The  de- 
gree of  accuracy  required  In  the  analysis  and  the  wing  plattform  shape  de- 
termines the  number  and  position  of  these  horseshoe  vortices. 

For  a given  wing  planform  geometry,  the  spam-  lift  distribution  Is 
determined  by  vhe  local  streamwise  angle  of  attack  distribution  and  the 
sectional  2-D  llft-curve-slope  distribution.  The  unknowns  In  the  problem  are 
the  values  of  running  lift,  p^,  at  points  along  the  wing  span.  To  compute 
the  relationship  between  the  local  angles  of  attack,  at  points  along  the 

wing  expressed  as  the  vector  and  the  values  of  running  lift  at  these 

same  points,  given  as  the  vector  the  following  linear  relationship  Is 
used. 

{»,}■  {B-D 
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where  is  referred  to  as  the  aerodynamic  influence  coefficient  (AIC)  | 

matrix. 

1 

To  compute  the  elements  ^^y  we  begin  by  subdividing  the  wing  into  a 
series  of  panels  of  spanwise  width  h^,  as  shown  in  Figure  B-1.  In  this 
figure,  the  horseshoe  vortices  are  bound  to  a point  at  the  wing  panel  quarter- 
chord  line  intersection  with  the  panel  midspan. 

For  a symmetrical  planform  at  zero  yaw  angle,  it  is  always  possible  to 
treat  any  spanwise  airload  distribution  as  the  sum  of  two  airload  distributions, 
one  which  is  symmetrical  with  respect  to  the  planform  centerline  and  the  other 
which  is  anti symmetrical  with  respect  to  the  planform  centerline.  This 
feature  considerably  simplifies  computation  of  the  A^j  elements  since  the 
specification  that  the  lift  distribution  is  either  symmetrical  or  antisym- 
metrical  allows  the  AIC  computations  to  be  restricted  to  one-half  of  the  total 
wing.  This  wing  half  is  taken  to  be  the  left  wing  of  the  aircraft. 

The  AIC  matrix  is  given  by  the  following  expression 


where 

a^  = 2-D  lift-curve-slope  at  wing  position  (i). 

The  elements  of  are  defined  in  terms  of  the  geometrical  position 

of  the  panel  quarter-chord  points.  These  expressions  are  somewhat  lengthy 
and  are  not  reproduced  here.  The  reader  is  referred  to  Reference  4 for  a 
detailed  discussion  of  the  derivation  of  the  AIC's. 
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APPENDIX  C 

TWIST  CAUSED  BY  CONTROL  DEFLECTION  PITCHING  MOMENT 

The  deflection  of  a control  surface  o wing  results  in  two  effects 
that  change  the  local  angle  of  attack  of  the  wing.  The  first  of  these 
effects  is  an  apparent  angle  of  attack  due  to  control  deflection,  written 
as 


where  6 is  the  control  surface  deflection.  The  second  effect  of  the  deflec- 
tion of  a control  surface  is  a twisting  of  the  wing  caused  by  a change  in 
the  pitching  moment  due  to  control  deflection.  It  is  to  this  latter  ef- 
fect that  attention  is  now  turned. 

When  an  aileron  is  deflected,  a streamwise  pitching  mopiCiH  results. 

This  pitching  moment  T is  given  by 

Ty  • <«> 

On  a flexible  wing,  this  pitching  moment  will  cause  deformation  that,  in 
turn,  will  lead  to  changes  in  the  spanwise  lift  distribution.  The  compu- 
tation of  this  lift  distribution  is  necessary  to  the  detenrii nation  of  the 
roll  control  effectiveness  of  an  aircraft. 

As  part  of  the  present  study,  the  matrix  relating  control  surface  de- 
flection to  the  streaixwise  angle  of  attack  was  determined  for  a laminated 
composite  beam  model.  For  this  model,  the  bending  moment  and  twisting  moment 
vectors,  computed  at  the  structural  reference  points  of  the  left  hand  wing 
panels  are  given  as  [4] 


{M}^  = -q[sinA][lQ]^[hc^c^j]{6> 

(C-3) 

{T}^  = q[cosA][lQ]^[hc^C|j^^]{6} 

(C-4) 
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utilizing  the  unit  load  method,  together  with  laminated  beam  theory,  the 
angle  of  attack  distribution  {a^}  caused  by  the  pitching  moment  given  in 
Equation  (C-2)  is  given  by  the  expression 


/ 

0 

0 

{“4 ' “ 

h 

/mU 

qh 

EI(1  - kg)  cosA 

rr 

EI(1  - kg)  COSA 

r / 
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♦ [t]  1 

0 

h 

0 

kh 

. (C-5) 
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L 

Vf  + 

6J(1  - kg)  COSA 

r>  j 

where 

[m]  = - 

[I^][sinA] 

(C-6) 

and 

[t]  = [Iq3Ccosa] 

(C-7) 

The  matrix  [1^]  is  defined  in  Equation  (A-17). 

Combining  Equations  (C-3)  and  (C-4)  with  Equation  (C-5)  yields  a rela- 
tionship of  the  form 

{-e}  = {«}  (C-S) 


where 


[Eij]  = [y 


h(sinA)(tanA-q) 
EI(1  - kg) 


h(cosA)(l  - ktanA) 
GJ(1  - kg) 


(C-9) 


The  matrix  defined  in  Equation  (C-9)  reduces  to  that  defined  in  Reference 
4 (p.  88)  for  a metallic  wing. 
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THE  STIFFNESS  EXPRESSIONS  FOR  A 
LAMINATED  COMPOSITE  BOX  BEAM 

The  wing  structural  model  employed  throughout  this  report  is  described 
in  detail  in  Reference  1 and  is  based  upon  the  assumption  that  the  bending 
stiffness  and  torsional  stiffness  of  the  wing  structure  are  due  entirely 
to  the  presence  of  relatively  thin,  laminated  composite  cover  sheets.  Al- 
though all  of  the  bending  and  torsional  stiffness  is  assumed  to  reside  en- 
tirely in  the  laminated  cover  sheets,  it  would  be  a simple  task  to  add 
algebraically  any  additional  stiffnesses  that  might  arise  from  the  presence 
of  such  elements  as  spar  caps  or  flexible  webs. 

To  formulate  the  expression  for  the  equivalent  bending  stiffness,  El, 
and  torsional  stiffness,  GJ,  for  this  type  of  box-beam  structure,  it  is 
necessary  to  make  an  assumption  about  the  deformation  behavior  of  the  wing 
box.  The  well-known  Euler-Bernoulli  hypothesis  that  the  strain  due  to 
bending  varies  linearly  from  the  neutral  surface  of  the  box-beam  provides 
this  assumption  about  the  displacement  behavior  of  the  wing.  In  addition, 
it  is  further  hypothesized  that  chordwise  sections  of  the  wing,  perpendicular 
to  a beam  reference  axis,  are  rigid.  This  means  that  the  deformation  of 
the  wing  box  is  a function  only  of  a spanwise  coordinate,  y.  The  material 
behavior  is  assumed  to  be  linear  elastic.  An  additional  feature  of  the 
idealized  beam  model  is  that  the  reference  surface  of  the  box  beam  is 
taken  at  the  geometrical  middle  surface  of  the  wing  box.  This  convention, 
commonly  used  in  laminated  plate  theory,  is  more  convenient  than  the  alter- 
native of  locating  a modulus  weighted  centroid  of  each  beam  cross-section 
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(a  modulus  weighted  centroidal  axis  will  always  lie  in  the  neutral  surface}. 
If  the  laminated  box  beam  is  not  both  elastically  and  geometrically 
symmetrical,  the  middle  surface  will  not  be  the  neutral  surface  for  bending. 

A Cartesian  coordinate  system  is  defined  in  Figure  D-1.  The  x-y  plane 
corresponds  to  the  geometrical  middle  surface  of  the  wing  box.  The  posi- 
tive x-axis  is  rearward  on  the  wing  while  the  y-axis  lies  along  the  swept 
wing  span  and  is  coincident  with  the  chordwise  centerline  of  the  wing  box. 


Figure  0-1  - Composite  material  lamina  principal 
axes,  1 and  2,  with  respect  to  the 
wing  reference  axes,  x and  y. 

For  a single  orthotropic  lamina,  the  relationship  between  membrane 
stresses  and  strains  occurring  in  the  beam  laminae  can  be  written,  in  matrix 
notation,  as: 


yy 


xy 


^22  ^26 
^26  ^66 


yy 


'xy 


(D-1) 


The  terms  Q. • are  functions  of  the  orthotropic  engineering  constants, 

* w 

Q. .,  and  the  angle  e defining  the  orientation  between  the  lamina  principal 

• w 

axes  and  the  box-beam  reference  axes  [5].  The  lamina  lies  in  a plane 
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^22  = b 
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1=1 
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(D-8) 


®22  ' ^ 


^22  ^ 


1=1 


(0-9) 


B33  = 2b 


1=1 


(D-10) 


and  b represents  the  chordwise  width  of  the  box  beam.  The  summation  in  the 
equations  extends  over  the  N layers  of  composite  material.  The  constants 
and  6^  are  defined  in  terms  of  lamina  coordinates  with  respect  to  the 
middle  surface.  In  terms  of  the  lamina  thickness  t^  and  lamina  lower  and 
upper  position  coordinates  z.  and  respectively,  b^  and  6.  are  defined 
as  follows 


/^i+1 


z^dz 


and 


= 3 {t^  + 3t^.Z;  + 3z.t^) 


■■■/ 


zdz 


"i 


(0-11) 


= ^j(Zi+i  + ^i)/2  (D-12) 

The  tena  b^  represents  the  area  moment  of  inertia  of  a strip  of  material, 
of  thickness  t^  and  unit  width,  about  the  middle  surface.  As  such,  b^  is 
always  a positive  number.  The  term  6^  represents  the  first  moment  of  the 
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ar^a  of  the  same  stHp  about  the  middle  surface.  This  term  «.  Is  positive 
if  the  lamina  area  centroid  lies  above  the  middle  surface  and  negative  if 
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